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Two antenna capacitance problems are solved in this paper. The motivating problem 
was to determine the capacitance of a biconical antenna of circular cross section immersed in 
a magneto-ionic medium. This is solved by reducing it to the secondary problem of deter- 
mining the capacitance of a biconical antenna of elliptic cross section in an isotropic medium. 
The latter is solved using the sphero-conal coordinate system. 



1. Introduction 

Satellite and rocket borne experiments carried out 
in the ionosphere utilize antennas for various 
purposes in ways which require an understanding 
of the effects of the ionospheric environment on the 
antenna impedance. In particular, many experi- 
ments are operated at frequencies and altitudes such 
that the anisotropic properties of the ionosphere 
come strongly into play. Examples occur in the 
reception of cosmic noise in the topside of the 
ionosphere at frequencies near 1 Mc/s [Walsh et aL, 
1963] and in the use of impedance probes to measure 
ambient electron densities [Herman, 1964]. In terms 
of the well known ionospheric parameters 

X= (plasma frequency /operating frequency)-, 
F=gyro frequency /operating frequency, and 
Z= collision frequency /operating frequency 

it turns out that over a large part of the XY 2 plane 
the capacitive part of the impedance is dominant. 
This capacitance can be thought of as made up 
approximately of parallel capacitances across the 
space charge sheath, i.e., from antenna surface to 
undisturbed plasma, and between the two parts of 
the antenna at different potentials. This latter 
capacitance has been computed in the literature by 
neglecting the modification of the embedding plasma 
due to the sheath [Herman, 1964; Kaiser 1962; and 
Balmain 1964]. Kaiser has considered a biconical 
antenna and Balmain a cylindrical one, while 
Herman's result is rough enough to require only a 
thin symmetrical antenna with no further shape 
details. 

To carry out his analysis of the conical antenna 
Kaiser [1962] employed an approximate mathe- 
matical procedure whose accuracy is difficult to 
assess. The purpose of this work is to obtain an 
exact solution for the same physically simplified 
problem; that of a conical antenna in a plasma, 



neglecting the sheath. This clears the way for evalu- 
ating the effects of the sheath by removing the effects 
of the approximate mathematics. The latter are 
negligible when the antenna axis is oriented along 
the earth's magnetic field II, but become significant 
if the antenna orientation approaches 90° to H. 
Both the present results evaluated for thin cones 
as well as Kaiser's and Balmain's [1964] results are 
all quite different from those of Herman. The pri- 
mary reason is that Herman [1964] incorrectly as- 
sumed that the value of the integral over the antenna 
surface of the electric field vector E is not different 
when the antenna is in free space from its value when 
the antenna is immersed in a plasma. 

Interest in biconical antennas of circular cross 
section stems in part from the fact that a thin 
biconical antenna approximates a thin cylindrical 
antenna and is easier to treat mathematically. In 
an isotropic medium the capacitance per unit length 
(along a generator) of an unterminated biconical 
antenna has the pleasant property of being constant 
and is given by 



c= 



In tan • 



In 



(!) 



(i) 



for small @. Here f3 is the semiangle of the cone and 
€ is the dielectric constant. The capacitance of a 
biconical antenna of circular cross section but finite 
length I measured along the generators is then simply 

C=cl 

if fringing is neglected. 

The plan of this paper is to first work out the 
analogous formulas for the case of a cone of elliptical 
cross section. This is done in section 2. Then in 
section 3 the reduction of the capacitance problem 
for a biconical antenna of circular cross section in 
anisotropic plasma to the problem solved in section 
2 is carried out. This is an application to the present 
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Figure 1. Elliptic cone in isotropic medium. 




Figure 2. The sphero-conal coordinate system. 



problem of a procedure given in Landau and Lifsbitz 
[I960]. 

Since use is made in section 2 of the rather little 
known sphero-conal coordinate system, a brief 
discussion of the form of the coordinate surfaces, 
and their analytic representation is included, largely 
following the work of Kraus and Levine [1961]. 

2. Capacitance of Elliptic Cone in Isotropic 
Medium 

Consider the boundary value problem 

d 2 V , d 2 V , d 2 F 







(2) 



dX 2 ' dX 2 ' dX 2 
and V= ± V on the surface of the elliptical cone 

X\ -\r\ 2 Xl = X 3 Xf (3) 

as shown in figure 1 ; X 2 , X 3 are positive and we assume 
X 2 >1. 

It is clear from the symmetry of the problem 
that the equipotentials are elliptic cones and the 
lines of force meridian lines. The natural coordinate 
system in this case is the sphero-conal system. 
Hobson [1931] and Kraus and Levine [1961] have 
given clear and concise treatments. The sphero- 
conal system uses concentric spheres and two sets, 
of elliptical cones as the three sets of mutually 
orthogonal surfaces. The concentric spheres are 
given by 

X 2 +X 2 +X 2 =r 2 . (4) 

The two sets of cones are shown in figure 2 and are 
given by 



■yi V2 V2 



x\ 



xi 



M 2 V- 

=0 



=0, 



(5) 



0<^ 2 </l 2 <M 2 </ 2 , 



Z>A>0. (6) 

Here I and h are arbitrary constants and for each 
value of n, system (5) represents an ellipse for 
constant X 3 with its major axis in the XiX 3 plane. 
Similarly for each value of v system (6) represents 
an ellipse for constant X x with its major axis in 
the X x Xz plane. The relation between the (X ly 
X 2 , X s ) and (r, /x, v) systems is given by 



^i — hi A 2 = — ,_. .. > 



hi 



x 3 



h^l 2 -h 2 



lyl 2 -h 2 



(7) 



The cone of our problem given by (3) can be made 
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to fit the notation of (5) if for convenience we take 
/x 2 = l and equate coefficients to obtain 



/r= — — and r — 



(8) 



The coordinates (r, p, v) are still not very con- 
venient for our purpose. Kraus and Levine [1961] 
have defined an equivalent set of polar coordinates 
(r, 6, ip) to replace (r, /x, v) with ranges of values 

0<d<7r, 0<^<2tt, r>0. 

Here 6 is the angle between the positive X 3 axis and 
the generator drawn through an end of the minor 
axis of the ellipse described by (5). Similarly, <p is 
the angle between the positive Xi axis and the gen- 
erator drawn through an end of the minor axis 
described by (6). Figure 2 may be used as an aid 
to visualize these surfaces. The following relations 
can be established 



cos <p=t> cos 6- 



V? 



where 



X 1 =r(l-Jt 2 cos 2 ^) 1/2 cos^, 
X 2 —r sin 6 sin <p, 
.Y3=rcos^(l-F 2 cos 2 ^) 1/2 , 



(9) 



k' 2 - 



lr 



tan 2 6 cos 2 S 



I 2 l + tan 2 ecos 2 



andJfc 2 =l-&' 2 . 



The semiangles of e and e' subtended by the foci 
of the ellipses of systems (5) and (6) at the origin 
respectively are given by 



tan e = 



h 
(l 2 -h 2 ) 1/2 cos 6 



and tan e = 



(l 2 -h 



2V/2 



h cos <p 



We summarize the salient geometrical features of 
the sphero-conal system. The surfaces r=const. 
are spheres centered at the origin. The 6= const, 
surfaces are elliptic cones belonging to system (5). 
Similarly, <p= const, surfaces are elliptic cones belong- 
ing to system (6). The equations of these surfaces 
can be found by eliminating \x, v in (5) and (6) with 
the help of (9). Two such intersecting surfaces are 
shown in figure 2. The line indicated as "common 
generator" is a line of intersection. This common 
generator traces out the elliptic cones of system (5) 
by rotating around the X 3 axis and when this gener- 
ator lies in the X 2 X 3 plane it makes an angle 9 with 
the X 3 axis. This same generator traces out the 
elliptic cone of system (6) by rotating around the 
Xi axis making an angle <p with the X x axis when it 
lies in the XiX 2 plane. Finally, when fc=l, the 
system degenerates into the familiar spherical system 
in the following fashion. As k— >1, the elliptic cones 
given by 6= const, become more and more circular. 
In the limit when fc=l, e=0 implying that the foci 
have coalesced and we have cones of revolution. 



As for the <p= const, surfaces, as &— >1, the major 

axis becomes longer and when fc=l, e f =- meaning 

that the foci have moved out to infinity, and the 
elliptical cone degenerates into a plane. 

In terms of these coordinates the cone of our 
problem is the surface d=6 where 

cos o =[l+A 3 /\ 2 )]- 1/2 . (10) 

The Laplace and gradient operator are given by 






1 



r 2 (k 2 sin 2 6+k' 2 sin 2 <p) 



[(i-F cos 2 ey» A ^ (x-v cos 2 ey> 2 | X 



+(i—k' 2 cos 2 ey* 



dp 

(l-k' 2 cos 2 <p) 1/2 £ 



(11) 



_ d, . (l-k 2 cos 2 6) 1 ' 2 d 

V r dr r(k 2 sin 2 6+k' 2 sin 2 <p) 1/2 bd 

(l-k' 2 cos 2 <pY< 2 d . . 

^ r(k 2 sin 2 6+k' 2 sin 2 v) 1 ' 2 d<p K ' 

Considering the Laplace equation we note that in 
the present problem there is no variation of V with 
respect to r and <p. Integrating we get 

V(6)=A f (l-k 2 cos 2 6)-" 2 de+B. (13) 

Jet 

The constants A and B can he determined from the 
boundary conditions 



V(e n )=v , V(w-e )=-v . 

Solving; for A and B we get 



/•' 



V(6)=V 



F 






(14) 



(15) 



Here F is the incomplete elliptic integral of the first 
kind given by 

F(6, k) = P (l-k 2 sin 2 <p)- i/2 <k>, k 2 <l. (16) 

The surface charge density is obtained with the help 
of the gradient operator (12) applied to the potential 
function (15). For this problem Z>=D • \t — a and 



1 eV 



•(f-M) 



(k 2 sin 2 e +k' 2 sin 2 < P )- 1/2 . (17) 
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It is seen that, in spite of the <p independence of 
V, a is <p dependent. The elemental area is given by 

,. rjk 2 sin 2 d +k' 2 sm 2 <pY< 2 
d0= (1-^cos^)^ - **• (18) 

Integrating D over < <p < 27r, < r < 1 and dividing 
by 2V , yields the capacitance c per unit length 
measured in the radial direction 



2K'e 



f 2 ' {1-k' 2 cos 2 <f>)- 1/2 <k> 



F 



(f» k ) 



(19) 



In the above result K' is the complete elliptic integral 
of the first kind of modulus k f . Here a cone of unit 
length is interpreted as that portion of the cone 
enclosed by a sphere of unit radius. For the co- 
ordinate system we are using, cones of this descrip- 
tion are suitable because integration in r will be 
particularly simple. 



3. Capacitance Problems in Anisotropic 
Media 

Assuming the static superposed magnetic field to 
be in the z direction, the dielectric tensor (of a mag- 
neto-ionic medium) has the form 



K= a 2 at • (20) 



<*1 


— a 2 


0" 


«2 


<*i 











<*3- 



In particular when collisions are neglected 

i X iXY A 1 y 

tti=l— vj _y 2 > ^2= , _y 2 ana a 3 =l— A, 

where X=-f> Y=— -> co p and co H are plasma and 



gyrofrequencies. 

The capacitance will be determined using a quasi- 
statics approximation in which electric field E is 
derivable from a gradient V (which is valid, strictly 
speaking, at zero frequency) , but using the dielectric 
tensor components evaluated at the frequency of 
interest. This procedure is justifiable, since capaci- 
tance is essentially a near field quantity, so long as 
the dimensions of the antenna are not large com- 
pared to the wavelengths of the radiation. This 
will be the case in most regions of the X— Y 2 plane 
shown in figure 3. However, for certain angles 6 
between the propagation vector and the superposed 
magnetic field the index of refraction for the ordinary 
mode becomes very large (infinite in the collision 
free case) in regions 7 and 8 and similarly for the 




Figure 3. The X-Y 2 plane. 

extraordinary mode in region 3. The associated 
wavelengths are then small and near zone or quasi- 
statics approximations become quite questionable. 
Furthermore, the radiation resistance in regions 3, 7, 
and 8 becomes very large [Weil and Walsh, 1964] 
and from this point of view one also sees that a quasi- 
statics procedure (which cannot yield radiation) 
cannot be expected to be valid. This will be dis- 
cussed further in section 4. 

The reduction of the capacitance determination 
for objects in a medium described by a Hermitian 

K to a related capacitance problem in an isotropic 
medium is demonstrated in Landau and Lifshitz 
[1960] beginning with the field due to a point charge. 
One has 



and 



V-D=e8(x)8(y)d(z) 



D=e *K-E. 



(21) 



Defining an electrostatic potential E=— VV and 
substituting in the above equations, we note that the 
mixed derivatives cancel on account of the Hermitian 
nature of the tensor and the equation for the potential 
V becomes 



dx 2 ^dy 2 ~ i ~ q bz 2 e^ 



,2+^ + <Z 2 ^2=~— *(*)%)*(*) 



(22) 



where 



The transformation 






x =x, y =y, z 



(23) 



reduces (22) to a conventional Poisson's equation 



b 2 V &V d 2 V 



5(x')8(y')8(z') (24) 



dx' 2 ' by' 2 ' ds" €o a l2 
for a point charge of strength e/2«i in vacuum 
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(primed space). Thus 



V= 



■M£)b' + »'+S'"- m 



The above concept is applicable in capacitance prob- 
lems since the surface charge densities must scale in 
the same way as individual charges. 

The capacitance C of an object in an anisotropic 
medium may now be determined as follows. First, 
one determines the capacitance C in vacuum of the 
related object obtained by applying transformation 
(23). If one assumes that changing from unprimed 
to the primed system and vice versa the potential of 
the body does not change then the total charge re- 
siding on the surface must change in the same manner 
as for a point charge. Based on this we have 



C=qa£'. 



(26) 



The determination of C is the secondary problem 
referred to in the abstract. 

The procedure employed by Kaiser [1962] was to 
transform the potential function obtained in the 
primed system back to the unprimed system ; then 
to obtain the total charge by integrating the normal 
component of D over the surface of the original 
antenna. Such a procedure has some drawbacks in 
general. First, the recovery of the potential function 
in terms of the unprimed system can be quite for- 
midable. Second, D and E will not be parallel in 
the unprimed system thus introducing further 
complications. Both these difficulties are avoided 
by using (26). 

We now revert to our specific problem. Consider 
a biconical antenna of semiangle whose axis makes 
an angle a with the magnetic field and lies in the yz 
plane as shown in figure 4. 

The problem is to solve 



bx 2 + dy 2 ' 



H d 



d 2 V 



0> 



(27) 



subject to the condition that V=± V on the surface 
of the cone given by 

x 2j r(— z sin a+y cos a) 2 =(z cos a+y sin a) 2 tan 2 p. 

(28) 

Applying the transformation (23) to (27) and (28), 
we get the following related problem in vacuum 



d 2 V , d 2 V , d 2 V 



172=0 



(29) 



dx' 2 ' dy' 2 ' i>z' 
and V= ± V on the elliptic cone 

z /2 +(l-sin 2 a sec 2 /% ,2 +r/ 2 (l-cos 2 a sec 2 f3)z' 2 

— 2(q sin a cos a sec 2 /3)y'z'=0. (30) 




Figure 4. Biconical antenna in magneto-ionic medium. 
z' 




*\l 



Figure 5. Elliptic cone of related problem given by {31), 



First, we must reduce (30) to the normal form (3) 
by the usual diagonalization process. Defining a 
new coordinate system £, 77, f rotated through an 
angle 5 from x', y' ', z' , as shown in figure 5, we have 
for the equation of the elliptic cone 



where 



^2, Xa 



^+X 2 r, 2 =X 3 f 2 



(1 + g 2 )— sec 2 j3 (sin 2 a-\-q 2 cos 2 a) 



(31) 



{ [(1 + g 2 ) -sec 2 (sin 2 q+g 2 cos 2 a)] 2 +4g 2 tan 2 <3} 1/2 
+ 2 

(32) 

X2 goes with the upper and X 3 goes with the lower 
sign. If X 2 >1, the elliptic cone (31) has the same 
orientation as that shown in figure 1 if we associate 
(£, V, f) with (x u x 2 , Xz). One can then determine 
O using (10). On the other hand if X 2 <1, we see 
from (8) and (9) that & 2 >1 and k' 2 is negative 
which is not permissible. This difficulty can be 
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overcome by rotating the £, rj coordinates through 
D0° but keeping f fixed. Replacing £ by — i\ x and 
77 by £1 in (31) we get 



{?+xjij?=x;r 2 



where 

A2 A2 

Now (33) is in the desired form. Using (10) 



(33) 



X2=— and Xs=^* 



cos 



«-[>+©]"" 



= (1+X 3 )" 



(34) 



At this point it may be pointed out that part of 
the mathematical approximation of Kaiser [1962] 
lies in using, instead of the conal equipotential 
boundaries which correctly appear in the related 
problem, an approximate set of elliptic cylinder 
equipotentials. In the present analysis the correct 
boundaries are maintained. Approximations are 
put in only at the end in the final formula for capaci- 
tance to facilitate numerical computations for thin 
cones. 

Before we apply the results of section 2 to this 
problem, we need an expression for the length of the 
elliptic cone. Let us assume that the radial length 
of the cone in figure 4 is unity. The coordinates 
{zo, Vo, 20) of the tip of a generator are 

Zo=sin P cos <p 

y = sin a cos /3+cos a sin P sin <p 

2 =cos a cos /3— sin a sin <p sin p (35) 

where (p is the azimuthal angle measured around the 
axis of symmetry z x (this <p is different from the <p 
of sphero-conal system) and varies from 0— 2tt. If 
(xo, y' , z'o) be the tip of a generator of the elliptic 
€one then its length is given by 



rS=(tf+»S , +s?) lfl 






using (23) 



q 2 — 1 ~| 1/2 

- 2 (cos a cos (3— sin a sin p sin <p) 2 

0<<p<2w. (36) 



We note that the generators of the elliptic cone are of 
different lengths except when a=0. For other 
angles of inclination, we will use an average value 
defined by 

( r 'o)=^ f Q * Kb)<h. (37) 

The above will lead to elliptic integrals. For small 
cone angles cos /3~1 and ignoring the terms con- 



taining sin P we get 

(r£)« ( sin 2 a- 



COS' 



a\ 1/2 t 



(38) 



In this small angle case (t'o) is just the length of the 
cone axis. 

The following special cases will now be considered. 

Case (1) : Axis of cone coinciding with the direction 
of the magnetic field. Here 



ol=0j X 2 =l, X 3 



: q 2 tan 2 p, 
6 =tan- 1 (gtanj8),A:=l,i / =0 (39) 

cos 2 A 1/2 



, / . 2/9 ,cos 2 /3Y 
r =(sin 2 /H — -J- J 



Using (19) for the capacitance per unit length of an 
elliptic cone in vacuum and the scale factor (26) 
we get 

^ Teoa^q 2 sin 2 ff+cos 2 p)W 

l'°Hi)]| 

Using trigonometrical relations we get 

tan (I)-?- 1 ! (cot 2 /3+<z 2 ) 1/2 +cot/3}- 

For this special case the exact result has also been 
obtained by Kaiser [1962]. For small values of 
P and q the above formula becomes 



c= 



7re ai 



In 



W 



(41) 



Case (2) : Axis of cone perpendicular to magnetic 
field. Here o=90°, X 2 =g 2 , X 3 =tan 2 P and (r£)« 1. 
Here we have two subcases q 2 ^>\ and q 2 <Cl. If 

tan § 

tan O = 



and 



with 



2g / « a a 1 g 



(42) 
(43) 



F=cos2/?+ !H^, k ^ sin ^_^iy 



For small P this becomes 

7re -V«i«3 



c=- 



Ln 



(i)' 



(44) 



If # 2 <1, the rotation of coordinates discussed 
following (32) is used. 

tan o =tan p 
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and 



with 



c= 



2K'e () a ] (j 



. tan 2 j8 



For small /3 this becomes 



'(HO 



(45) 



1+^ 



7reoV«i«3 
As against these Kaiser [1962] obtains 

7T€oV«l«3 



(46) 



"■yy 



(47) 



Case (3): Axis of cone making intermediate 
angles with magnetic field. 



2,vV/2 



<r )«^sin 2 oH — - 2— 1 



so that 



2/T 



600:3 



8*(J-M) 



/ . 2 , cos 2 a \ 1/2 .... 

^sin-a+— ^-J • (48) 



6 — 




J L 



J L 



Figure 6. 
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Capacitance as a function of inclination a. 
X"=0.835, F 2 =3.33, /3=1°. 
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Figure 7. Capacitance as a function of inclination a. 

A r =0.784, ys=0.19fi, /3 = 1°. 
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Figure 8. Capacitance as a function of inclination a. 

.Y=3.33, F2=o.835, = 1°. 



4. Numerical Results and Discussion 

In practical situations where only cones of small 
angle are of interest the following approximations 
are valid 

¥^\ and & /2 ->0 
so that 



K' < 



and 



2p(!-* 0> *)«|lntan(Jf) 



The resulting simplified capacitance formula 

,1/2 



7re a 3 



8 In cot Q) 



( . 2 cos 2 aY 



(49) 



will be utilized in numerical calculations. It must 
be noted that a further approximation of setting 
6 = /3 is not valid. This follows by considering (39) 
and (42) . In the former 6 = tan -1 (q tan 0) and in the 

latter ^ =tan~M 1. Since these values depend 

on q, one cannot set 6 ~I3. 

As we propose to compare present results with 
those of Kaiser's [1962], his formula is given below 
for convenience 



/ . 2 .cos'aV" 

7r€ a 3 l sin a-\ ^— 1 



if 4(^ 2 sin 2 a+cos 2 a) "^ 

q \qp[l + (cos 2 a+q 2 sin 2 a) 1 ' 2 } J 



(50) 



In evaluating numerical results, it is necessary to 
know the ranges of values of q 2 for which formula 
(49) is valid. It may be recalled that we were re- 
quired to solve the following partial differential 
equation 



b 2 V b 2 V 2 d 2 F 
dx^dy 2 ^^ dz 2 



(27) 



with Diriclilet conditions. In a classical sense we 
have a well posed problem provided (27) is elliptic, 
i.e., 5 2 >0. In the primed space we have simply to 
solve Laplace's equation with V given on real 
boundaries. To take into account complex values 
of q 2 corresponding to nonnegligible electron-neutral 
collisions in the plasma one can analytically con- 
tinue the above expressions for c, valid for real q 2 , 
into the complex q 2 plane, avoiding the necessary 



branch cuts. However, the physical meaning of the 
results becomes doubtful as one nears the negative 
real axis where q 2 <^0 and the equation in fact be- 
comes hyperbolic. This is because q 2 <^0 corresponds 
to the collision free case for X and I 72 combinations 
which lie within regions 3 and 7 of figure 3. In these 
regions, as previously pointed out, because of the 
short wavelengths in certain propagation directions, 
near zone (i.e., quasi-static) approximations are not 
Aalid. In addition, antenna impedance depends 
intimately on the current distributions and these are 
to be expected to be highly modified in regions 3 and 
7 from those of an electrically short antenna. 

For these reasons we would therefore not attribute 
physical significance to the result obtained by letting 
q 2 be negative, namely that there is an imaginary 
part in c corresponding to an apparent radiation 
resistance and real loss of energy from the antenna 
to the collisionless plasma in regions 3 and 7, in addi- 
tion to the radiation loss which is found by using 
the full Maxwell equations. 

In accord with the above remarks we will consider 
only regions of # 2 >0 for numerical results. Repre- 
sentative values are shown graphically in figures 
6 to 8 corresponding to regions 6, 1, and 5, respec- 
tively. In region 5, X>1, consequently capacitance 
is negative meaning that the antenna will be induc- 
tive. From the computed results we see that 
Kaiser's approximation is in agreement with the 
present results only for angles of inclination a 
around zero and deviates for increasing angles of 
inclination, reaching a 20 percent difference in some 
cases. 
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